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Abstract 

We report on the first experimental verification of the Zurek-Kibble scenario in an isolated su- 
perconducting ring over a wide parameter range. The probability of creating a single flux quantum 
spontaneously during the fast normal-superconducting phase transition of a wide Nb loop clearly 
follows an allometric dependence on the quenching time tq, as one would expect if the transition 
took place as fast as causality permits. However, the observed Zurek-Kibble scaling exponent 
a = 0.62 ± 0.15 is two times larger than anticipated for large loops. Assuming Gaussian winding 
number densities we show that this doubling is well-founded for small annuli. 

PACS numbers: 11.27. +d, 05.70.Fh, ll.10.Wx, 67.40.Vs 
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The Zurek-Kibble (ZK) scenario [H |2j [3] proposes that continuous phase transitions take 
effect as fast as possible i.e. the domain structure after the quenching of the system initially 
reflects the causal horizons. This proposal can be tested directly for transitions whose 
domain boundaries carry visible topological charge. In this letter we shall show that, with 
qualifications, this scenario is strongly corroborated by the behaviour of superconducting 
loops, for which the topological defect is a fluxoid i.e. a supercurrent vortex carrying one 
magnetic flux quantum $ = h/(2e) ~ 2.08 x lCT 15 Wb. 

The basic scenario is very simple. Consider a planar low-T c superconductor in which 
a hole has been made of circumference C. In the Meissner state the order parameter for 
the superconductor is a complex field ip with phase cf>, ip = pe 1 ^, where \p\ 2 measures the 
density of Cooper pairs. On quenching the system from the normal to superconducting 
phase, causality prevents the system from adopting a uniform phase. If, on completion of 
the quench, we follow the periodic phase <f>(x) (mod 2ir) along the boundary of the hole 
(co-ordinate x), we can define a winding number density: n(x) = d(f>(x) / dx / (2ir) . The total 
normalized magnetic flux through the hole is, in units of $o> the winding number: 



where A<ft is the change in <p. In the absence of an external magnetic field, on average 
(n) = 0, but it will have non-zero variance (An) 2 = (n 2 ), which is what can be measured in 
terms of the probabilities f± m to trap ±m flux quanta: (n 2 ) = J2m=-oo m2 fm- According 
to Ref.[2], on completion of a thermal quench having a given inverse quench rate tq = 
—T c /(dT/dt)T=T c , the phase is correlated over distances 27r£, where £ was predicted to 
depend alio metrically [2j on the quench time tq: 



£, also called the ZK causal length, is defined in terms of the cold correlation length £o and 
the Ginzburg-Landau relaxation time r of the long wavelength modes. The ZK scaling 
exponent a is determined by the static critical exponents of the system and, in the mean- 
field approximation, a = 1/4 PQ. If we make the further assumption that there is a random 
walk in phase on a scale 2ttC, then, for a hole of radius r, circumference C ^> 2-7r£, 




(1) 




(2) 
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(n 2 >^- = f (VT (3 ) 

For small rings with C < 27r£ the likelihood of seeing two or more units of flux is small 
and (n 2 ) « f +1 + /_j = / 1( the probability of single fluxoid trapping. It is plausible to 
extrapolate Eq.(|3]) to: 

/i«<n 2 >«f f^l ■ (I) 



showing allometric behaviour of /i with the same exponent. We note that the ZK argument 
makes no assumptions about the rest of the superconductor, equally valid for the phase 
change along the inner circumference of an annulus as it is for the phase change around a 
single hole in a superconducting sheet. In 2003 the first experiment with superconducting 
loops [1] was performed to test Eq.(|3]). The experiment consisted of taking an isolated array 
of thin-film wide rings and making it undergo a forced phase transition by heating it above 
its superconducting critical temperature and letting it to cool passively back towards the 
LHe temperature. Once the thermal cycle is over, the rings are inspected by a scanning 
SQUID and the number and polarity of any trapped fluxoids determined. These rings, of 
amorphous MosSi thin films, had thickness almost one order of magnitude smaller than the 
low temperature London penetration depth. Although this provides favorable conditions for 
thermally activated phenomena it drastically increases the likelihood that nucleated vortices 
escape through the ring walls during the fast quench. In fact, the experimental outcome 
was totally at variance with the allometric scaling above. However, the prediction Eq.pJ 
presupposes that we can ignore the contribution to the flux from the freezing in of thermal 
fluctuations of the magnetic field [S] and the results of [I] could be explained in terms of the 
freezing of thermally activated fluxoids in a similar spirit to Ref.j5]. 

In this paper we shall present results from a new experiment with high-quality Nb film 
rings with r = 30 /im, two times thicker than their low temperature London penetration 
length Xi^Nb'i the film thickness and composition were chosen to reduce the thermal activation 
of fluxoids and, at the same time, 'the washing out' of fluxoids generated by the conventional 
causality mechanism. In our case the contribution A/i to the probability of finding a unit 
of flux from thermal magnetic field fluctuations is approximately [5]: 
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A/i < {k B T c )r^ 2 « 6 x lO" 4 , (5) 

and can be safely ignored. We therefore look for scaling behaviour in tq. 

Here a different way of counting both the number and the polarity of generated defects 
has been adopted. It is based on the detection of the persistent currents J s circulating 
around a hole in a superconducting film, when one or more fluxoids are trapped inside the 
hole. The circulating currents screening the bulk of the superconductor from the trapped 
flux induce a magnetic field H in the volume around the ring, such that J s = VxH. By 
placing a Josephson tunnel junction ( JTJ) along the perimeter of the hole in the area where 
this field passes, any trapped fluxoid will result in a modulation of the JTJ critical current, 
similar to the effect of an external field applied perpendicular to the ring. Indeed, this 
method is strongly inspired by the results found investigating the effects of a transverse field 
on Josephson junctions of various geometries [6]. The geometry of our experiment is sketched 
in FigjTJ the black wide ring is a 200 nm thick Nb film, which also acts as the common base 
electrode for two JTJs whose top electrodes are depicted in grey. The JTJs have the shape of 
gapped annuli and the bias current is supplied in their middle point: to our knowledge, this 
geometrical configuration has never been realized before and is characterized by a peculiar 
magnetic diffraction pattern: for small magnetic fields the critical current increases both for 
positive and negative field values. The original purpose of having two counter electrodes 
on the base ring was that any screening current circulating on the outer ring circumference 
will preferentially affect the outermost JTJ, and vice versa for the screening current on 
the inside of the ring. Since the persistent currents due to trapped flux mainly flow in 
the inner ring circumference [7j, in this experiment we only used the innermost JTJ. The 
layout shown in FigjTJ with a few key differences, bears remarkable topological similarity to 
the one used in a series of experiments by us to demonstrate the ZK scaling behaviour of 
Eq. ([3]) in annular JTJs [8j [9j [TOj [TTJ [12] . The most obvious difference in the design is the 
inclusion of the two junction counter electrodes on top of the ring-shaped base electrode. 
The second change is the removal of a small section of the full annular junction to leave a 
gapped annular junction with the purpose of avoiding fluxons created inside the JTJ at the 
Josephson phase transition. Should any be produced, they will simply migrate through the 
junction extremities driven by the applied bias current needed to overcome eventual pinning 
potentials. This leaves the experiment only sensitive to the fluxoids produced in the ring at 
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FIG. 1: Sketch of a superconducting loop (black) used as a base electrode for two gapped 
Nb/AlOx/Nb annular Josephson tunnel junctions (whose top electrodes are in gray). The ring 
inner and outer radii are r = 30 and R = 50/nn, respectively, while the top electrodes width is 
5/um. 

the phase transition. 

As with our previous experiments, the present one relies on a fast heating system, obtained 
by integrating a Mo resistive meander line on the 4.2mmx3mmx 0.35mm Si chip containing 
the ring with the Nb/AlOx/Nb JTJs. The quench time tq could be continuously varied over 
more that four orders of magnitude (from 20 s down to 1 ms) by varying the width and the 
amplitude of the voltage pulse across the integrated resistive element. In order to determine 
the quench time with high accuracy, the ring temperature was monitored exploiting the 
well known temperature dependence of the gap voltage of high-quality Nb/AlOx/Nb JTJs 
already described in Ref.|10j. After each ring thermal quench the critical current of the 
innermost JTJ is automatically stored and an algorithm has been developed for the counting 
of the trapped fluxoids. Finally, all the measurements have been carried out in a magnetic 
and electromagnetically shielded environment. During the thermal quenches all electrical 
connections to the heater as well as to the JTJs were disconnected. While more details on 
the measurement setup and on the fabrication process can be found in Ref . [12] and Ref . [2] , 
respectively, an extensive description of the chip layout, the experimental setup and the 
system calibration will be given elsewhere |1 5] . 

The experimental results shown in Fig. [2] were obtained using a ring with inner and 
outer radii, r = 30 /im and R = 50 /im, respectively. Similar samples have shown the 
same behaviour. We note that wider rings prevent fluxoids from tunneling out of the ring, 
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although their smaller normal self-inductance L n makes fluxoid formation energetically more 
unlikely. In our case, the field energy E = §l/2L n associated with a single flux quantum 
$o is several orders of magnitude larger than the thermal energy ksT c /2 at T„|lb'j. 

Magnetostatic numerical simulations implemented in the COMSOL Multiphysics 3D Elec- 
tromagnetics module showed that, when a single flux quantum is trapped in such rings, the 
radial magnetic field induced by the circulating currents at the ring inner border is as large 
as lA/m, a value easily detectable by the JTJ. For our rings, the number of trapped flux- 
oids was small, usually no more than one; indeed we measured the probability of trapping 
a single up-fluxoid (field up) f + \ and a single down-fluxoid 

Fig. [2]shows on a log-log plot the measured frequency f\ = f + i + /_i = (n + i ±n_i)//V = 
rii/N of single fluxoid trapping, obtained by quenching the sample ./V times for each value of 
a given quenching time tq, ri\ being the number of times that one defect or one anti-defect 
was spontaneously produced. N ranged between 250 and 300 and n\ was never smaller 
then 10, except for the rightmost point for which ri\ = 5. The vertical error bars gives the 
statistical error f\/\/n\. The relative error bars in tq amounting to ±10% are as large as 
the dot's width. As expected we had n + \ m n_i, but slightly larger than n_i indicating the 
presence of a small residual stray field in our apparatus (see inset of Fig. 2). 
To test Eq.(2), we have fitted the data with an allometric function fi = aTQ b , with a and 
b as free fitting parameters. An instrumentally weighted least-mean-square fit of fi vs. tq, 
represented by the continuous line in Fig. [2j yields a = 0.18 ± 0.02 (taking tq in ms) and 
b = 0.62 ±0.15. The large fit correlation coefficient R 2 = 0.987 indicates that the allometric 
behaviour is reliably confirmed, however the scaling exponent b is about two times larger 
than expected for large loops. 

A doubling of the large-loop ZK exponent for small loops has a possible explanation in the 
framework of the Gaussian correlation model introduced in Ref . [12] in which it was assumed 
that the winding number n(x) is a Gaussian variable until the transition is complete, whereby 
all correlation functions are determined by the two-point correlation function g(xi — X2, C) = 
(n(xi)n(x2)). As a result [T7] : 

(n 2 ) = / / (n (xi) n (x 2 )) dxidx 2 = 2C / g(x,C)dx. (6) 
Jo Jo Jo 

For C 27r£, we can assume a correlation function of the form g(x, C) = 

g(x/2rt^C/2nO/(2rt0 2 , so that: 




FIG. 2: Log-log plot of the measured frequency f\ of trapping single fluxoid versus the quenching 
time tq for a Nb ring having inner radius r = 30 /im, outer radius R = 50 /im and thickness 
d = 200 nm ~ 2A^ ) Arfc. Each point corresponds to hundreds of thermal cycles. The vertical error 
bars gives the statistical error, while the relative error bars in tq amounting to ±10% are as large 
as the dots' width. The solid line is the best fit to an allometric relationship f\ = CLTq h which 
yields a = 0.18 ± 0.02 (taking tq in ms) and b = 0.62 ± 0.15. For comparison purposes, the dashed 
line is the prediction of Eq.di}])) with x (see text) set to 0.7 to fit in ordinate scale. 



fi « 2—- g(x, C/2n0dx « 2 ^—-) g(0, C/2^). (7) 



provided g(x, C/2ir^) is analytic at x = 0. This suggests that Eq.(|4j) should be replaced by 
the scaling behaviour: 

M£)'-«(5)'Gr- 

with a proportionality constant k of the order of unity. To buttress this suggestion, it is not 
difficult to show that, in the Gaussian approximation, the value of (n 2 ) along a small ring 
in a 2D superconductor is proportional to the area enclosed by the ring |15j . 
This doubling of the scaling exponent in Eq.Q has the price of coming with a lower prob- 
ability, but leaves us with some freedom with the ZK prefactor. Indeed, the value of the 
prefactor a obtained from the allometric best fit of the experimental data in Fig|2] is about 
k = 4-5 times larger than the predicted value (r /£,o) 2 y/r = 0.04 obtained using the values 
r = 30 /im, £o ~ 30 nm and tq = 7ih/16kBT c « 0.16 ps and taking tq in ms. As a bound we 
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only expect agreement in the overall normalization of the prefactor a to somewhat better 
than an order of magnitude, largely confirmed by experiment. We point out that the de- 
pendence of the prefactor a on the ring width remains to be investigated both theoretically 
and experimentally [We note that, if Eq.(|3]) were true, then (n 2 ) would be pa 0.6, i.e., 20 
times larger than the experimental value for tq = O(10ms), say] 

In the opposite case of large circumferences, C ^> 27r£, g(x,C) is controlled by the 
correlation length £ of the winding number at the time of unfreezing and does not depend 
on C, i.e., the effect of periodicity for large rings is small. With g{x) = g(x/27r£)/(27r£) 2 on 
dimensional grounds, we justify the random walk assumption of Eq.(|3]), 

(n 2 ) = -j / g(z)dz pa -= / g(z)dz = -=, (9) 
< Jo < J Q 2vr£ 

with x — 0(1). The dashed line in Figj2]is the prediction in Eq.Q with x set to 0.7 to fit 
to the ordinate scale. 

In summary, our experiment shows reliable scaling behaviour of the form Eq.(j3]) for the 
creation of a single fluxoid, with scaling exponent 0.62 ± 0.15. This is obviously at variance 
with the extrapolation Q of the Zurek prediction Eq.([3]) to small rings, for which we expect 
a = 0.25. We have suggested that it be given by Eq.(j8]) with twice the exponent, as a 
consequence of the Gaussianity of the Cooper pair field phase (before truncation by back- 
reaction), an assumption supported in a slightly different context by the behaviour of JTJs 
in an external field [12J. As Gaussianity permits the instabilities from which defects form to 
grow as fast as possible, in general it gives the same scaling exponents as the ZK scenario 
for large systems. With this qualification we see our result as providing strong support for 
Zurek-Kibble scaling over a wide range of quenching time tq. We stress that this experiment 
is the only one to date to have confirmed the Zurek-Kibble causality scenario for single 
isolated superconducting rings (as distinct from Josephson junctions). Further experiments 
to investigate the transition to the random walk regime and the effect of the ring width are 
planned. For example, a test of Gaussianity is that fx < 0.5 for all values of C [T2] . 
Given that the original ZK scenario was posed to demonstrate the similarity in the role of 
causality at transitions in the early Universe and in condensed matter systems we have seen 
that the fmiteness of the latter systems requires careful disentangling from the underlying 
principles before we can draw any quantitative conclusions. 

In the same vein we conclude with a speculation concerning the small-annulus JTJ exper- 
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iments [9J [TUl HH 12] that, hitherto, have been the only superconductor experiments to show 
scaling behaviour. In that case also, the observed exponent was twice that anticipated from 
long annuli [8j. However, in the case of JTJs there is an ambiguity in their fabrication that 
is sufficient to double the exponent, according as the 'proximity effect' enables otherwise 
subcritical behaviour of the Josephson current density to dominate near the transition \X0\ . 
We had assumed that this was the reason for the discrepancy. We shall now reexamine these 
earlier experiments with the above analysis in mind. 

The authors thank P. Dmitriev for the sample fabrication and testing, A. Gordeeva for 
useful discussions and M. Aaroe for the help at the initial stage of the experiment. 
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